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We present 350 /im photometry of all 17 galaxy candidates in the Lockman 

^ : Hole detected in a 1.1 mm Bolocam survey. Several of the galaxies were previously 
detected at 850 /xm, at 1.2 mm, in the infrared by Spitzer, and in the radio. 

^ ! Nine of the Bolocam galaxy candidates were detected at 350 /im and two new 

CN ! candidates were serendipitously detected at 350 fim (bringing the total in the 

(y-^ I literature detected in this way to three). Five of the galaxies have published 

^ , spectroscopic redshifts, enabling investigation of the implied temperature ranges 

\^ I and a comparison of photometric redshift techniques. 

^ ' A = 350 /im lies near the spectral energy distribution peak for z ~ 2.5 ther- 

o 



mally emitting galaxies. Thus, luminosities can be measured without extrapo- 
lating to the peak from detection wavelengths of A > 850 /im. Characteristically, 
^ I the galaxy luminosities lie in the range 1.0 — 1.2 x 10^^ L0, with dust temper- 

atures in the range of 40 K to 70 K, depending on the choice of spectral index 
and wavelength of unit optical depth. The implied dust masses are 3 - 5 x 10^ 
Mq. We find that the far-infrared to radio relation for star-forming ULIRGs 



C/3 



H ' systematically overpredicts the radio luminosities and overestimates redshifts on 

the order of Az ~ 1, whereas redshifts based on either on submillimeter data 
alone or the 1.6 /im stellar bump and PAH features are more accurate. 
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Subject headings: galaxies: high-redshift — galaxies: starburst — submillimeter 

1. Introduction 

Surveys at submillimeter and millimeter wavelengths have detected hundreds of galaxy 
candidates by their thermal dust emission. The galaxies (hereafter referred to as submil- 
limeter galaxies) characteristically have redshifts z > 1 and inferred luminosities of L ~ 10^^ 
L0 and star formation rates of 10^ Mq per year (assuming dust heating by young stars). 
Such enormous luminosities and star formation rates, or accretion rates in the case of super- 
massive black hole growth, make submillimeter galaxies strong candidates for the progenitors 
of massive galaxies at the current epoch. Clearly, it is crucial to characterize the spectral 
energy distributions (SEDs) where their emission peaks (A = 350 fim for 40 K dust at a 
redshift of z = 2.5), measure their redshifts and luminosity functions, determine their power 
sources, and integrate them into theories of galaxy formation. 

Although submillimeter galaxy SEDs peak at a few hundred microns for all but the 
highest redshifts, A > 850 fim surveys have been most successful at detecting galaxies be- 
cause of the lower atmospheric noise and greater transmission, and less stringent telescope 
surface requirements. Most of the detections have been low signal-to-noise ratio (just over 
thresholds of 3-4 cr), necessitating multiwavelength confirmation. Furthermore, the SEDs 
have been extrapolated shortward from 850-1200 //m over the peak, or between 850-1200 iim 
and the far-infrared, to estimate dust temperatures, luminosities, and star formation rates. 
Clearly, 350 /im photometry can confirm galaxy candidates and sample the SEDs near their 
peaks for more precise inferences of physical parameters. Similarly, because of the difficulty 
in obtaining spectroscopic redshifts of large numbers of highly obscured galaxies, various 
photometric redshift estimation techniques have emerged, notably based on the far-infrared 
to radio luminosity relation in ULIRGs (CariUi & Yun 1999; Yun & Carilli 2002) and the stel- 
lar continuum bump in the infrared for Spitzer-detected galaxies (Egami et al. 2004; Sawicki 
2002). However, despite the difficulty, candidate spectroscopic redshifts have been obtained 
for ~ 73 galaxies (Chapman et al. 2005). Thus, with well-determined dust-emission SEDs, 
including 350 /im, photometric techniques can be compared to spectroscopic redshifts. 

In this paper, we present 350 /im photometry of all 17 submillimeter galaxy candidates 
from the Bolocam Lockman Hole survey (Laurent et al. 2005), some with previous 850 /im 
and 1200 /im detections. Two new galaxy candidates were serendipitously detected at 350 
/im, bringing the total number of 350 /im-discovcrcd galaxies to three (Khan et al. 2005). 
We combine these data with infrared and radio data to derive improved luminosities, ex- 
plore the range of implied dust temperatures and spectral indices, and compare photometric 
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redshift techniques. This comparison is timely with the imminent launch of the Herschel 
Space Observatory (scheduled for August 2007), which will detect thousands of galaxies at 
far-infrared and submillimeter wavelengths, but for which spectroscopic redshifts will be at- 
tainable for only a small fraction. Throughout the paper, a cosmology of Hq — 70 km s~^ 
Mpc~^, flM = 0.3, and fl\ = 0.7 is assumed. 

2. The 350 //m SHARC II Galaxy Survey 

Observations at multiple submillimeter wavelengths are vital both to confirm the Bolo- 
cam sources (as 6 false detections are expected from Monte-Carlo simulations) and to make 
photometric redshift and temperature estimates. The 350 fim. SHARC-11 observations com- 
bined with the Bolocam 1.1 mm galaxy survey provides a flux density ratio that is strongly 
dependent on redshift for a given temperature. This is because the rest wavelength cor- 
responding to the observed wavelength of 350 //m with SHARC-II is near the peak of the 
grey-body spectrum (for a z ~ 2 galaxy at 40 K), and Bolocam's 1.1 mm observations climb 
the steep (/3 fa 1.5) modified Rayleigh- Jeans side of the spectral energy distribution. 

Follow-up observations of each of the Bolocam Lockman Hole galaxy candidates (Lau- 
rent et al. 2005) were taken with the Submillimeter High Angular Resolution Camera 
(SHARC II) at the Caltech SubmiUimeter Observatory. The observations were taken on 
three observing runs: 2004 March-April, 2005 January, and 2005 February. The brightest 
Bolocam sources (1 and 2) were observed over 8 hours of total integration time during the 
2004 March- April run, although most of the run was lost due to poor weather. Bolocam 
sources 1 and 16 were observed over 18 hours of integration time during the 2005 January 
run, again with much of the run lost due to poor weather. The 2005 February run was 
characterized by much better weather, and all the Bolocam sources except for 1, 5, 8 and 16 
were observed over 35 hours of integration time^. When combined with the observations of 
Bolocam sources 5 and 8 by Kovacs et al. (2005), the entire Bolocam sample was observed 
over these observing runs. 

Observations with SHARC 11 were taken in the point source observing mode, with a 
Lissajous (parametric sinusoidal curve) scan pattern using the SWEEP command of tele- 



^The primary weather measurement correlated with the SHARC II mapping speed is the CSO 225 GHz 
heterodyne, narrowband, "tipper tau" monitor, which measures the zenith atmospheric attenuation. The 
2004 March- April, 2005 January, and 2005 February Lockman Hole observations yielded T225GHz ranges and 
75th percentiles of T225GHz = 0.046 - 0.087, T75% = 0.076, T225GHz = 0.044 - 0.120, T75% = 0.093 and 
r225GHz = 0.030 - 0.074, t75% = 0.047, respectively. 
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scope. The Lissajous pattern was scanned in altitude and azimuth, with amphtudes of 30" 
and 20", respectively. When combined with the 2.6' x 1.0' SHARC II field-of-vicw and 9" 
FWHM instrument beam size, this resulted in a uniform coverage region of 95" x 18", with 
a border of additional coverage (60" x 40") outside this region. Each observation had a fixed 
length of 10 minutes to ensure uniform coverage even on individual scans. Integration times 
and the resulting depths of each of the SHARC II fields are listed in Table 1. 

The reduction of the raw SHARC II data was accomplished with the use of the " deep" 
cleaning utility of the Comprehensive Reduction Utility for SHARC II (CRUSH^). Obser- 
vations of pointlike galaxies, quasars, protostcllar sources, Hn regions, and evolved stars 
were used to construct pointing models for each of the observing runs. Observations of the 
pointing sources were taken with a scan strategy identical to that of the science fields. A 
subset of the pointing sources were used for flux density calibration, with reference 350 /im 
flux densities obtained from the SHARC II website^. 

Source extraction was performed on the CRUSH-cleaned maps, with each map (corre- 
sponding to a single Bolocam candidate) consisting of all of the individual scans co-added 
together. The algorithm was begun by doing a cut on the uniform coverage region, defined 
as the set of pixels for which the coverage is > 60% of the maximum per-pixel coverage. The 
uniform coverage region is a contiguous region in the center of each map. Next, an RMS 
in sensitivity units (the flux density of each pixel times the square root of the integration 
time for that pixel in units of mJy s^^^) was computed in the uniform coverage region. This 
RMS is valid for the entire uniform coverage region since variations in coverage have been 

1/2 

accounted for by the t/ coverage normalization, where ti is the total integration time for 
pixel i. All pixels with coverage- normalized flux densities exceeding 3 a ("hot pixels") were 
flagged as potential sources. Then hot pixels were grouped into multi-pixel sources by making 
the maximal group of adjacent hot pixels, including those within y/2 pixels (i.e., diagonally 
adjacent). The peak flux density, right ascension and declination of the source candidates 
were computed by centroiding two-dimensional Gaussians on the groups. The uncertainty 
in the flux density of each source is given by the pixel-to-pixel RMS at the centroid location 
of the source. 



^http:/ /www.submm.caltech.cdu/^sharc/crush/index.htm 
^http: / /www. submm.caltech.edu/ ~sharc/ 
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Table 1. SHARC II Photometry and New Galaxy Candidates 
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Note. — SHARC II detections and 3 cr upper limits at each of the Bolocam sources, in order of descending brightness at 1.1 mm. 
The Bolocam sources in parentheses correspond to SHARC II detections well outside of the Bolocam 2 cr positional error circle, 
and are therefore believed not to be associated with the Bolocam source. Two SHARC II detections from Kovacs et al. (2005) are 
also included. The SHARC II source numbers are listed in order of 350 /im S/N. 
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3. Positional Uncertainties 

The large beam sizes of submillimeter and millimeter wave instruments (31", 14", 11", 
and 9" FWHM for Bolocam, SCUBA, MAMBO, and SHARC II, respectively) makes it 
difficult to identify likely optical and radio counterparts to the galaxy candidates. Despite the 
large beam sizes, however, individual sources can be centroided to much higher precision than 
the quoted beam size. To help constrain this issue of source matching between the various 
surveys, a positional error circle was estimated for each of the submillimeter and millimeter 
band detections. For the Bolocam detections, Monte Carlo simulations were performed 
by injecting sources into the timestream and running the reduction pipeline and source 
extraction algorithm. This simulation was repeated for a range of source flux densities. The 
resulting centroiding error as a function of flux density (5.4 - 9.1") was added in quadrature 
with the RMS telescope pointing error (9.1") to yield a range of 2 a positional error circles of 
21 - 26". A similar approach was used to estimate the centroiding errors for SCUBA (Scott 
et al. 2002) and MAMBO (Greve et al. 2004), yielding 2.2 - 10.4" and 1.2 - 4.3" respectively. 
When added in quadrature to the quoted pointing errors (4" and 3", respectively), this yields 
2 a positional errors of 9.2 - 22" and 6.5 - 13", respectively. As the centroiding error as a 
function of flux density for the SHARC II observations was not available, values of 3.0 - 4.0" 
were empirically determined from the SHARC II centroiding flts of the Bolocam sources, 
which, when added in quadrature with the pointing error of 3.8", yields 2 a positional errors 
of 9.8 - 11". These error circles were used to correlate the sources between the different 
surveys to find coincident detections. 



4. Results 

4.1. SHARC II 350 /xm Detections 

Postage stamp images of each of the SHARC II fields are shown in Figure 1. Each 
image has been cropped to 60" x 60", centered on the Bolocam source positions (dotted 
circle). The SHARC II source candidate fist is presented in Table 1, where the sources are 
listed in order of Bolocam source number. Seven Bolocam galaxy candidates were detected 
by SHARC II at > 3a (Bolocam 1, 2, 3, 6, 14, 16, 17). Two of these sources (Bolocam 
3, 6) were found to have two SHARC II counterparts. Two Bolocam candidates (5 and 8) 
were observed by Kovacs et al. (2005) and are also included in the list. 3 a upper limits 
are given for each of the Bolocam fields with no positive detections. Note that the fiux 
density uncertainties in the last column of Table 1 include uncertainties of ~ 20% due to 
calibration error (as determined by the dispersion of the calibration source flux densities). 
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The correlation between the two source lists based upon the positional error circles and 
detection offsets is lower than expected, as only 2 of the 7 Bolocam sources with a single 
SHARC II counterpart have errors within 1 a (5 are expected from a normal distribution). 
This may be due to underestimating the Bolocam pointing error (cf. § 3 and Laurent et al. 
2005). 

An additional two sources were detected in the survey (in the fields of Bolocam 9 and 
16); they are not associated with the Bolocam sources because their locations are well outside 
of the Bolocam positional error circles (see § 3). While 0.8 false detections are expected from 
Gaussian statistics in the SHARC II survey (given the 60% uniform coverage cut and the 9" 
SHARC II beam), these serendipitous sources nevertheless may be real. SHARC II sources 3 
and 4 were detected at high significance at 350 /im and lie in regions of positive flux density 
in the Bolocam map, with S/N ratios of 2.9 and 2.1, respectively. 

4.2. SHARC II / Bolocam Correspondence 

In addition to the Bolocam 1.1 mm and SHARC II 350 fim detections, existing mul- 
tiwavelength coverage (submillimeter, radio, infrared, optical, and X-ray) of the Lockman 
Hole was used to identify likely counterparts and characterize coincident sources. A detailed 
description of each of these surveys is found in Appendix A. A comprehensive summary 
of the counterparts to the Bolocam detections (including the coverage by each survey) is 
listed in Table 3. Given the large size of the Bolocam beam, identifying likely counterparts 
requires a certain amount of judgement. Detailed maps of the sources can be found in Figure 
2. Additional notes on individual objects are discussed below. 

Bolocam.LE. 1100.1 - We conclude that Bolocam source 1 is hkely to be a submillimeter 
galaxy given the coincident Bolocam, MAMBO, and SHARC II detections. In addition, 
highly plausible radio, Spitzer, and faint optical counterparts exist. 

Bolocam.LE.1100.2 - The SHARC II detection falls 16" to the ENE, but contains both 
radio sources in its error circle. The two 20 cm radio sources (Yun et al. 2005; Biggs & 
Ivison 2006) are unrelated due to their separation. We treat each of the radio sources 
separately when fitting photometric redshifts. The coincident northern radio source with 
bright optical counterparts has a low photometric redshift (see next section) and is likely to 
be a low-redshift galaxy as the SDSS survey concludes. Bolocam and SHARC II may also 
be detecting the Biggs & Ivison (2006) southern radio source, which has a very faint optical 
counterpart. 

Bolocam.LE. 1100.3 - The position of this 6.0 mJy Bolocam source is at the edge of 
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Fig. 1. — Postage stamp images (in sensitivity units of mJy s^/^) of each SHARC II field, 
centered on the Bolocam source positions. Images have been cropped to 60" x 60". The 
dotted circles represent the 2 a positional error circles of each Bolocam galaxy candidate. 
Overlayed are contours of signal-to-noise, with levels of 1, 2, 3, 4, 5, and 6. Contours of 
positive SHARC II detections (S/N > 3) are thicker than the others for clarity. Bolocam 
sources 5 and 8 are not included, as they were observed by Kovacs et al. (2005). 
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Fig. 2. — Plot of coincident multiwavelength detections of all 17 Bolocam sources, centered 
on the Bolocam positions. Sizes of the Bolocam, SCUBA, MAMBO, and SHARC II circles 
correspond to the 2 a positional error circles. Black circles correspond to VLA radio de- 
tections by Yun et al. (2005), Ivison et al. (2002), Ciliegi et al. (2003), and Biggs & Ivison 
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the good coverage region of the SCUBA survey. All of the optical counterparts to the 
radio sources are relatively bright (22 - 24 magnitude), although curiously, the SDSS catalog 
classifies the optical counterpart to the northeast radio source as a star. (Given the radio and 
350 fim SHARC II counterparts, we conclude that the SDSS classification may be incorrect.) 
Given that both of the SHARC II sources have radio, Spitzer, and optical counterparts, each 
(or both) are likely candidates as submillimeter galaxies and each could contribute to the 
flux density of the Bolocam source. Note that when estimating photometric redshifts where 
source confusion may be present (for this Bolocam source and elsewhere), no attempt was 
made to partition the Bolocam flux density among multiple submillimeter sources (due to 
the large uncertainties in position). 

Bolocam.LE. 1100.5 - Four radio detections fall within the Bolocam positional error 
circle: one within the SCUBA and SHARC II error circles and just outside the edge of the 
MAMBO error circle, one on the edge of the MAMBO error circle, and the other two near 
the edge of the Bolocam positional error circle. The southwest radio source has 5-color SDSS 
photometry and is classified as a galaxy (extended). The fact that three of the radio sources 
lie outside both the SCUBA and SHARC II error circles makes them unlikely to be the 
correct counterpart of the submillimeter detections. We therefore choose the northeast radio 
source to be the more likely counterpart (which is confirmed by the fact that Chapman et al. 
2005 were able to obtain a spectroscopic reshift for this submillimeter galaxy at this radio 
position, as discussed in Appendix B). 

A MAMBO detection located just outside the Bolocam error circle was not detected by 
SHARC II or SCUBA (3 a upper limit of 35.4 mJy), but has a Chapman et al. (2005) 
spectroscopic redshift of 1.956. Due to the large size of the Bolocam beam, the Bolocam flux 
densities and positions may be influenced by source confusion. 

Bolocam.LE. 1100.6 - Two SHARC II counterparts fall within the Bolocam positional 
error circle, each with radio counterparts (with the eastern source containing two radio 
counterparts). The SDSS survey classifles both the radio source associated with the west- 
ern SHARC II source and the eastern radio source associated with the eastern SHARC II 
source as galaxies. Each of the three radio sources may be contributing to the Bolocam flux 
density due to source confusion. We treat each of the radio sources separately when fltting 
photometric redshifts. 

Bolocam.LE. 1100.8 - Given the fact that the southern radio source (of the pair of 
two radio sources oriented N-S) has Bolocam, SCUBA, MAMBO, SHARC II, and Spitzer 
detections, along with a spectroscopic redshift, we conclude that this Bolocam source is real. 
Nevertheless, the northern radio source cannot be ruled out as a galaxy also contributing to 
the submillimeter fluxes. 
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Bolocam.LE.1100.14 - This 4.4 mJy Bolocam detection is likely influenced by source 
confusion, given three closely spaced submillimeter sources (SCUBA sources 1, 4, and 8, with 
the latter two lying near, but outside of the Bolocam positional error circle). The location 
of the southern radio source relative to the Bolocam position is greater than the 16" radius 
used for the Spitzer counterparts catalog, and thus no Spitzer data is available. This set of 
coincident sources is the most likely counterpart to the Bolocam source. Just outside of the 
Bolocam error circle hes another 850 //m SCUBA detection (LE850.4) to the northeast, with 
a coincident Ivison et al. (2002) 20 cm radio source as well as pubhshed 3.6, 4.5, 5.8, and 
8.0 /im Spitzer counterparts. 

The SCUBA source that coincides with Bolocam source 14 (LH850.1) is also the bright- 
est SCUBA source and has been the target of many pubhshed multi-wavelength studies. In 
addition to the extensive radio, infrared, optical and X-ray surveys discussed in Appendix 
A, a faint (K ~ 23.5) galaxy counterpart was positively identified (Lutz et al. 2001) at the 
radio position. The source was found to be extended (20-30 kpc), clumpy (on subarcsecond 
scales) and very red (/ - K > 6.2). 

Bolocam.LE.1100.16 - We conclude that this Bolocam source is in fact a submillimeter 
galaxy, given the large number of multiwavelength detections and a radio source with a 
confirmed spectroscopic redshift. 

Bolocam.LE.1100.17 - This 4.0 mJy Bolocam detection is hkely infiuenced by source 
confusion, given two nearby submillimeter sources. The 850 iim SCUBA, 1.1 mm MAMBO, 
and 850 //m SHARC II coincident detections to the northeast are the hkely counterparts 
to Bolocam source 14 and are discussed in detail in Bolocam.LE.1100.14. The SDSS cat- 
alog curiously classifies the northwest radio source as a star. We conclude that both the 
southeast and northwest radio sources associated with the Bolocam source may be submil- 
limeter galaxies, given the large number of multiwavelength detections. While a confirmed 
spectroscopic redshift exists near the southeast radio source, self-consistent photometric red- 
shifts and multiple optical counterparts at the radio position suggest that the spectroscopic 
redshift may not correspond to the radio / submiUimeter sources (see § 5.4). 

4.3. SHARC II Non-Detections 

From extensive Monte-Carlo simulations of the Bolocam data set (Laurent et al. 2005), 
6 false detections (Poisson distributed) are expected. This represents a large fraction (6/17) 
of the overall source catalog and is a consequence of the relatively low 3 a detection threshold 
used in the source detection algorithm. Eight of the Bolocam sources (4, 7, 9, 10, 11, 12, 
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IS, and 15) were found to show no secure counterparts at 350 /im, although two of the 
sources (Bolocam 9, 12) exhibit flux densities just below the 3 a detection threshold (with 
a coincident radio detection for Bolocam source 9). Here we describe each of the SHARC II 
non-detections of the Bolocam sources in detail. 

Bolocam.LE. 1100.4 - A single radio counterpart (Biggs Sz Ivison 2006) lies near the edge 
of the MAMBO positional error circle, with an SDSS classification of the optical counterpart 
as a galaxy. While well within the Bolocam positional error circle, the location of the radio 
source relative to the Bolocam position is greater than the 16" radius used for the Spitzer 
counterparts catalog, and thus no Spitzer data is available. Given the coincident Bolocam 
and Mambo detections, along with a lack of SHARC II and SCUBA detections, this source 
could possibly be a very high redshift galaxy (z > A), such that the SED falls below the 3 a 
detection threshold of the SCUBA 850 ^m survey. 

Bolocam.LE. 1100.7 - The lack of multiwavelength observations makes it difficult to 
determine whether this Bolocam detection is real (and associated with the coincident radio 
detection). The lack of Spitzer and SHARC II counterparts to the radio source, however, 
leads us to believe that the Bolocam source may be a spurious detection. 

Bolocam.LE. 1100.9 - It is interesting to point out that the SHARC II upper limit in the 
Bolocam error circle is just below the 3 a detection flux density threshold and coincides with 
the radio position. The lack of multiwavelength observations makes it difficult to determine 
whether this Bolocam detection is real. The presence of Spitzer counterparts and a possible 
dim SHARC II detection, however, leads us to beheve that the Bolocam source may be real. 

Bolocam.LE. 1100. 10 - Given the lack of counterparts, we conclude that there is little 
evidence to suggest that this detection represents a submillimeter galaxy and is likely a 
spurious detection. 

Bolocam.LE. 1100. 11 - The lack of more multiwavelength data makes it difficult to 
determine whether this Bolocam detection is real. While lacking a SHARC II detection, the 
radio source with an SDSS classification as a galaxy (extended object) leads us to beheve 
that the Bolocam source may be real. 

Bolocam.LE. 1100. 12 - A portion of the Bolocam error circle hes outside the deep Ivison 
et al. (2005) optical Subaru R-band field. Similar to Bolocam source 9, we point out that the 
SHARC II upper limit in the Bolocam error circle is just below the 3 a detection fiux density 
threshold. The lack of more multiwavelength data makes it difficult to confirm whether this 
Bolocam detection is real. 

Bolocam.LE. 1100. 13 - This Bolocam source lacks 850 //m SCUBA detections (although 
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a SCUBA source is located just outside of the Bolocam positional error circle). We suggest 
that there is little evidence that the Bolocam detection represents a submillimeter galaxy 
and is likely a spurious detection. 

Bolocam.LE.1100.15 - The lack of MAMBO and SHARC II counterparts makes it diffi- 
cult to confirm the Bolocam detection. Nevertheless, the SDSS classification of the position 
coincident with the radio source as a galaxy (extended) leads us to believe that the Bolocam 
source may be real. 

4.4. Submillimeter Spectral Energy Distributions 

The submillimeter spectral energy distributions (SEDs) of the coincident SHARC / 
Bolocam detections is shown in Figure 3. Five of the 17 Bolocam galaxy candidates (5, 8, 
14, 16, 17) have spectroscopic redshifts from Chapman et al. (2005, see Appendix B). In 
order to properly compare the SEDs, it is necessary to shift each of the SEDs to a common 
redshift. Thus, each observed SED was brought to a redshift of 2.0 (the mean redshift of the 
five Bolocam galaxies) using the spectroscopic redshifts. The composite SED of these five 
Bolocam galaxies can be seen in Figure 4. 

In addition to redshifting the SEDs to align their rest wavelengths, a cosmological 
dimming term was applied by assuming a flat {Qk = 0), Q\ = 0.7 cosmology. Finally, to 
account for variations in their intrinsic brightnesses, we normalize the flux densities of these 
five Bolocam galaxies by tying together their SEDs at the observed Bolocam wavelength of 
1.1 mm. To account for the spread of the redshifted wavelengths of the 1.1 mm Bolocam 
observations, the fiux densities were normahzed to the Laurent et al. (2005) model based 
on the observations cited in the Blain et al. (2002) paper. The model assumes a single dust 
temperature of 40 K (/? = 1.6) and is ovcrplotted as a solid line in Figure 4. Note that only 
the Bolocam observations are constrained to pass through this model. 

Upon inspection, we find that at least four of the five Bolocam galaxies with spectro- 
scopic redshifts (5, 8, 14, 16) exhibit very similar SEDs in the submillimeter and millimeter 
wavelengths. They are modeled adequately by the 40 K composite SED based on nearby 
dusty IRAS galaxies, high redshift submillimeter galaxies, gravitationally lensed high-redshift 
galaxies, and high redshift AGN. The Bolocam galaxy 17, however, appears to peak at a 
much higher wavelength and lower flux density than the others. We believe that there is 
enough source confusion to question whether the Chapman et al. (2005) redshift for this 
galaxy (z — 0.689) corresponds to the SED shown (see § 5.4). If the spectroscopic redshift is 
valid, the SED is modeled much better by a T = 20 K (/3 = 1.0) grey-body dust spectrum. 
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Fig. 3. — Submillimeter SEDs for all Bolocam galaxies detected in at least one other sub- 
millimeter waveband. The solid line represents a fit to a modified blackbody spectrum using 
the parameters of Laurent et al. (2005, See § 5). 
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Fig. 4. — Composite submillimeter SED for the five Bolocam galaxies with spectroscopic 
redshifts (Bolocam 5, 8, 14, 16, and 17). The observed SEDs have all been redshifted to 2; = 

2.0 using the spectroscopic redshifts, corrected for cosmological dimming, and their Bolocam 

1.1 mm fluxes normalized to the Laurent et al. (2005) model. The dotted points represent 
Bolocam source 17 assuming a redshift of 2.4 (the median photometric redshifts from § 5.2). 
The solid line represents the Laurent et al. (2005) model (T = 40 K, /3 = 1.6) based on the 
composite SED of Blain et al. (2002, and references therein). The dashed line represents the 
same model, with T = 20 K and /3 = 1.0. 
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5. Redshifts 
5.1. Introduction 

With the multiwavelength photometry of the Bolocam sources, we fit photometric red- 
shifts using various models based on different portions of the SED. Photometric redshifts 
based on the far-IR-to-radio correlation were derived using the models of Carilli & Yun 
(1999) and Yun & Carilli (2002). The shape of the submillimeter and millimeter part of the 

spectrum was also fit without the radio points, assuming a blackbody emission spectrum 
modified by a dust cmissivity term (Wiklind 2003; Laurent ct al. 2005). A brief description 
of each of the models, fitting methods, and the redshift results are discussed in the next 
section. 



5.2. Redshift Techniques 

This section attempts to briefiy describe each of the five photometric redshift techniques 
used in this paper and the results of the fits when applied to the Bolocam galaxy candidates. 
The following section will compare the relative merits of each of the photometric redshift 
fitting techniques and discuss the results of redshift distributions. 

1) FIR-Radio Spectral Index - Carilli & Yun (1999) used the semianalytic, linear rela- 
tionships derived by Condon (1992) between the massive star formation rate and the radio 
synchrotron luminosity and far-IR dust emission from active star-forming galaxies to show 
that the spectral index between these two frequencies, af^^, is a well behaved function of 
redshift: 

° = -0.24 - [0.42 X (a^adio - ttsubmm) X log(l + z)], (1) 

where we adopt the standard value in Condon (1992) of -0.8 for ctradio, and a value of +3.2 
for asubmm (au average of the spectral indices between 270 and 850 GHz for M82 and Arp 
220. The relation is believed to be a result of relativistic electrons accelerated in supernova 
remnants (producing synchrotron radiation) and dust heated by the interstellar radiation 
field (with a thermal peak of ~ 380//m for a galaxy with z — 2 and T = 40 K). Photometric 
redshifts determined using only the Bolocam and radio flux densities are listed in Table 2. 
Redshift results from Bolocam sources with multiple radio counterparts are listed using the 
higher S/N detection in the case of coincident detections by independent surveys or are listed 
together in the case of multiple counterparts detected by a single group. 

The error bars listed in Table 2 (and elsewhere throughout this paper) were obtained 
from Monte-Carlo simulations of the flts and represent statistical errors due to measurement 
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uncertainty in the flux densities. The flux densities at each observed wavelength were varied 
about their mean value assuming a Gaussian distribution of flux errors. Each Monte-Carlo 
SED was then flt to the photometric redshift models with a standard, least-squares min- 
imization fitting routine. Each simulation was repeated 1000 times, with the error bars 
quoted being the minimum-length 1 a confidence intervals from the resulting histogram of 
redshifts. It should be noted that these confidence intervals represent only the statistical 
goodness of flt and that uncertainties in the templates themselves are expected to dominate 
the photometric redshift errors. 

2) Entire FIR-Radio SED - Yun & Carilli (2002) utilized the entire Far-IR to radio 
spectral energy distribution to estimate photometric redshifts and SFRs. The redshift tem- 
plate is based upon the theoretical models of thermal dust emission, thermal bremsstrahlung 
(free-free) emission, and nonthermal synchrotron emission for dusty starburst galaxies. Pho- 
tometric redshift flts of the flve Bolocam galaxy candidates (5, 8, 14, 16, 17) with spectro- 
scopic redshifts (Chapman et al. 2005) are shown in Figure 5, with best-flt redshifts (and 
errors) also listed in Table 2. The solid lines in the flgure represent the best flt spectrum to 
the submillimeter, millimeter, and radio point shown. The dotted line represents a second 
flt using the Yun & Carilli (2002) model, this time flxing the spectroscopic redshift and 
normalizing (varying only the SFR) to the submillimeter points. 

3) Modifled Blackbody - Wiklind (2003) found that observations of local ULIRGs exhibit 
a remarkably low dispersion in the far-IR to millimeter wavelengths (A > SO/um), independent 
of whether the power source of the thermal emission is due to AGN or intense star formation. 
Wiklind (2003) flt a simple blackbody emission spectrum (modifled by a dust emissivity term) 
to sample of 37 local ULIRGs from Klaas et al. (2001): 



where B,^{T) is the Planck function evaluated at dust temperature, T, and frequency, u, and 
Ti, is the optical depth of the dust: 



Wiklind (2003) made no assumption about the Wien side of the spectrum, as only the 
submillimeter (> 450 fim) and millimeter points were flt. Using the best-flt parameters 
from Wiklind (2003): (3 = 1.8, uq = 1.2 x lO^^ (250 fim), and = 68 K, we flt 
photometric redshifts to the SHARC II 350 /xm, SCUBA 450 and 850 /xm, Bolocam 1.1 mm, 
and MAMBO 1.2 mm flux densities of the galaxies detected in our Bolocam survey. The two 
parameter flt (redshift and overall flux density normalization) yields redshifts for 9 of the 17 
bolocam galaxies with > 2 submillimeter/millimeter points. Seven of the Bolocam galaxies 



oc e^B^{T) oc [1 - exp {-t^)]B^{T), 



(2) 




-19- 



100 
10 

1 

0.1 
0.01 

0.001 
0.0001 
100 
10 

1 

0.1 

0.01 

0.001 
0.0001 
100 
10 

1 

0.1 

0.01 

0.001 
0.0001 




Bolocam 5 



Bolocam 8 



Bolocam 14 



Bolocam 16 




3olocam 17 



Yun & Carilli (2002) Photo- 




Bolocam 



SFR 

Msol/yr 



SFR 

Msol/yr 



5 


O AO.S 
^.^0.2 


660™ 


2.611 


718°: 


8 


4 7" 

4: . f 0.4 


758™ 


3.036 


871S 


14 


O 1 0.2 
^. J- 0.2 


726™ 


2.148 


661"= 


16 


O qO.2 
^.^0.2 


54711 


2.142 


624a: 


17 


O QO.3 
<'.0o.3 


445:; 


0.689 


158^; 



10° IQi 102 103 10^ 105 
Wavelength, /^m 



10° 101 102 



103 10* 105 



Fig. 5. — Photometric redshift fits to five of the Bolocam galaxy candidates, using the model 
of Yun & Carilli (2002). The fits were limited to the sub millimeter, millimeter and radio 
points shown. The solid line represents a two parameter fit, allowing both z and SFR to 
vary. The dashed line is the result of fixing the redshifts to the spectroscopic redshifts of 
Chapman et al. (2005), normalized to the submillimeter points. The best fit redshifts and 
SFRs, along with their error bars (from Monte Carlo simulations) are listed for each source. 
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(Bolocam 7, 9, 10, 11, 12, 13, 15) have no counterpart in the submiUimeter/miUimeter and 
one (Bolocam 4) has detections only at 1.1 and 1.2 mm, which is an insufficient wavelength 
spread in order to properly constrain the galaxy redshift using this two parameter model. 
Redshift results for each of the 9 Bolocam galaxies are listed in Table 2. The best fit models 
to the five galaxies with Chapman et al. (2005) redshifts (Bolocam 5, 8, 14, 16, 17) are shown 
in Figure 6. 

Similar to the method of Wiklind (2003), Laurent ct al. (2005) created a composite SED 
of nearby dusty IRAS galaxies, high-rcdshift submillimcter galaxies, gravitationally lensed 
high-rcdshift galaxies, and high-rcdshift AGNs (Blain ct al. 2002, and references therein), 
and found fit parameters of T = 40 K, z/q = 3700 GHz, and /3 = 1.6 for Equation 2. Redshift 
results for each of the 9 Bolocam galaxies are listed in Table 2. The best fit models to the 
five galaxies with Chapman et al. (2005) redshifts (Bolocam 5, 8, 14, 16, 17) are shown in 
Figure 6. 

4) Near IR Stellar Bump - Egami et al. (2004) used the multiband imaging capabilities 
of the IRAC and MIPS IR cameras of the Spitzer Space Telescope to observe 38 VLA radio 
sources in the Lockman Hole. They classified the resulting IR SEDs into two types: those 
showing a clear near-IR stellar continuum hump at a rest wavelength of 1.6 /im (due to the 
minimum opacity of the H~ ion at 1.6 /xm from photo-detachment and free- free transitions, 
which results in a local maximum in the the SEDs of cool stars, Sawicki 2002), and those 
with a featureless power-law continuum (from AGN). Wc fit the Spitzer IR, counterparts of 
the Bolocam galaxies with each of these spectra, using the Egami ct al. (2004) models of a 
cool ULIRG Arp 220 (from Silva et al. 1998) and a warm (dominated by an AGN) ULIRG 
Mrk 231. Only two fit parameters were used: the redshift, and an overall normalization. 
Photometric redshifts were fit for each of the 12 Bolocam galaxies (Bolocam 1, 3, 5, 8, 9, 10, 
11, 13, 14, 15, 16, 17) with > 2 IR Spitzer points. Three Bolocam galaxies (Bolocam 2, 6, 
12) were outside of the field surveyed by Egami et al. (2004). Two galaxies (Bolocam 4, 7) 
have Spitzer counterparts, but due to the high density of Spitzer sources in the field, they 
could not be uniquely associated with the Bolocam sources (because of a lack of another 
coincident detection in the submillimcter and/or radio). The best-fit photometric redshifts 
for both Arp 220 and Mrk 231 (fitting only the Spitzer near- and mid- infrared points) are 
shown in Figure 7, with the resulting redshifts also listed in Table 2. 

5.3. Compcirison of Photometric Redshift Techniques 

Comparing the results of each of the photometric redshift techniques with the spec- 
troscopic redshifts of Chapman et al. (2005) yields widely varying degrees of agreement. 
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Fig. 6. — Photometric redshift fits to five of the Bolocam galaxy candidates, using the 
model SEDs of Laurent et al. (2005, solid line) and Wiklind (2003, dotted line). The fits 
were limited to only the submillimeter and millimeter points shown. The best fit redshifts, 
along with their error bars (from Monte Carlo simulations) are listed for each source. 
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Fig. 7. — Photometric redshift fits to five of tlie Bolocam galaxy candidates, fitting the 
Spitzer near- and mid-infrared points to Egami et al. (2004) models of a cool ULIRG Arp 
220 (soUd line, from Silva et al. 1998) and a warm (dominated by an AGN) ULIRG Mrk 231 
(dotted line). The (sub)millimeter and radio points were not used in the fit and are shown 
only for reference. The dashed line represents a hnear combination of the two models). The 
best fit redshifts along with their error bars (from Monte Carlo simulations) are listed for 
each source. 
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Histograms of redshift errors for each of the photometric redshift models (other than CariUi 
& Yun 1999, in which coincident radio detections were treated separately) are shown in Fig- 
ure 8. The histogram from fitting models of Arp 220 and Mrk 231 to the Spitzer IRAC and 
MIPS observations are the fits that result in the lowest residual (Arp 220 for Bolocam 5, 
8, 14, and 16, and Mrk 231 for Bolocam 17). 

The Yun & Carilli (2002) model (as well as CariUi &; Yun 1999) yields systematically 
high photometric redshifts compared to the spectroscopic redshifts by Chapman et al. (2005). 
The comparison of model SEDs for photometric (solid line) and spectroscopic (dotted line) 
redshifts in Figure 5 suggests that SED data points at the extreme ranges of wavelength 
coverage strongly influence the model fit and that the systematic tendency to derive a high 
redshift is primarily driven by the lower than expected radio continuum flux density. This 
is supported by the fact that flts to only the submillimeter and millimeter-wave points yield 
much more accurate photometric redshifts (see below). This is perhaps not surprising, given 
recent evidence (Chapman et al. 2005) suggesting a large degree of dispersion in the radio- 
to-far-lR correlation at higher redshift. Nevertheless, it is unlikely that the Yun & Carilli 
(2002) model template, which is derived from the ensemble average of 23 infrared luminous 
galaxies in the local universe, is systematically biased by radio bright objects because >98% 
of all FIR-selected galaxies follow the well known and tight radio-FIR correlation, indepen- 
dent of FIR luminosity (Yun, Reddy, & Condon 2001). Aside from Bolocam source 17 whose 
spectroscopic redshift by Chapman et al. (2005) appears suspect (see § 5.4), these compar- 
isons suggest that the observed radio continuum in Bolocam galaxies is 2-5 times fainter 
than predicted by the synchrotron flux densities (which dominate thermal brehmsstrahlung 
by a factor of ~ 13 at 20 cm) from the low-redshift ULIRGs for which the local FIR-radio 
correlation was derived (see further discussions in 6.2). 

In contrast to the far-IR-to-radio correlation photometric redshift techniques, both the 
Wikhnd (2003) and Laurent et al. (2005) modifled blackbody curves correctly estimate the 
redshifts for three of the five Bolocam galaxies with spectroscopic redshifts (within the 
uncertainties of the photometric redshift techniques). The strength of the submillimeter / 
millimeter-only photometric redshift technique is twofold. First, while abandoning the radio 
points limits the number of points (as few as two, in some cases) to which we may fit a 
model, we ensure that the physics that dominates the region of the spectrum to which we 
are fitting is directly relevant to star formation-heated dust emission. Synchrotron radiation, 
by contrast, is dependent upon high energy electrons streaming through interstellar magnetic 
flelds, whose properties may vary as a function of environment (e.g. inverse- Compton losses 
for cosmic rays with higher CMB energy densities at high redshift - see § 6.2). Indeed, these 
galaxies are typically at least an order of magnitude more luminous than the low-redshift, 
infrared-luminous galaxies from which the FIR-radio correlation was derived. 
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Second, having the Bolocam 1.1 mm flux densities on the Ray leigh- Jeans side of the 
spectrum and the SHARC II 350 /xm flux densities near the peak of the SED makes the 
350 /xm/l.l mm flux density ratio a strong function of redshift. This can be seen in Figure 
9, which shows the flux density ratios between various wavebands based on the Laurent 
et al. (2005) model SED. The importance of the SHARC II 350 fim waveband is apparent. 
For intermediate to high {z < 5) redshifts, the SHARC II flux density drops rapidly as 
a power law (~ due to the hotter components of dust) with redshift on the Wien 

side of the spectrum, while the millimeter-wave climbs up the steep Ray leigh- Jeans portion 
of the SED. Photometric redshifts using the 450 or 850 iim wavebands are less sensitive 
than the 350 /im / 1.1 mm wavebands. Extending this analysis shows the discriminatory 
power of the BLAST and Herschel space telescope 250 fim. bands in conjunction with a 
millimeter waveband, although the far-IR waveband begins to probe a range of hotter dust 
temperatures. (The flux density ratio of these wavebands may not be well correlated, as 
discussed in § 6.1.) BLAST (Balloon-borne Large- Aperture Submillimeter Telescope, Devlin 
et al. 2001) is a balloon-based instrument which incorporates a 2-meter primary mirror and 
is equipped with large-format bolometer cameras operating at 250, 350, and 500 /im which, 
when complete, will provide the first sensitive large-area (0.5-40 deg^) submillimeter surveys 
at these wavelengths. The bolometer arrays are prototypes of the Spectral and Photometric 
Imaging Receiver (SPIRE) focal plane cameras for the Herschel satellite (Griffin, Swinyard 
& Vigroux 2001), which will further investigate the formation and evolution of AGNs and 
star formation in high redshift submillimeter galaxies. 

It is important to note that while both the Wiklind (2003) and Laurent et al. (2005) 
modified blackbody curves both produce reasonably accurate photometric redshifts (and 
produce nearly identical SEDs on the Rayleigh- Jeans portion of the spectrum), they make 
substantially different assumptions about the dust properties: Tjust and f3 are 68 K and 1.8 
for the Wiklind (2003) model SED and 40 K and 1.6 for the Laurent et al. (2005) model 
SED. This points to the degeneracy of the dust temperature and the grain emissivity index. 
While the shapes of the submillimeter SEDs are reasonably modeled by either dust model 
and thus predict photometric redshifts with some accuracy, essentially no information about 
the dust temperatures can be inferred. In fact, representing the dust SED with two (or more) 
components produces similar values (and thus similar redshifts), with both temperatures 
lower than that of the single dust temperature model (Wiklind 2003) . 

The photometric redshifts determined by fitting the Spitzer infrared observations yield 
redshifts that are equivalent to both the Laurent et al. (2005) and Wikhnd (2003) model 
SEDs. This confirms the conclusions of Egami et al. (2004), in which starburst-dominated 
galaxies (" cold" ) show remarkably similar SEDs in the infrared. The resulting photometric 
redshifts are highly sensitive to the 1.6 /im continuum hump (and PAH features), with a 
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sharp minimum in . (The ARP 220 fits may be biased towards particularly good fits as 4 
of the 5 galaxies with spectroscopic redshifts lie between 2 < 2; < 3, which is optimal for the 
8 [xm PAH feature to be shifted into the observed 20 /xm IRAC waveband.) AGN-dominated 
("warm") galaxies also show very similar SEDs, but lacking a strong continuum feature in 
the infrared, are subject to larger redshift fitting uncertainties; a brighter, higher redshift 
galaxy is characterized by a similar shape in the infrared portion of the SED as a cooler, low- 
redshift galaxy, with the Wien side of the spectrum well-modeled with a power-law (Blain 
et al. 1999). 



5.4. Bolocam Source 17: Spectroscopic Misidentification? 

We bring special attention to photometric redshift analysis of Bolocam source 17 (cor- 
responding to SCUBA source 8), as the Laurent et al. (2005), Wiklind (2003), and Spitzer 
IR models are consistent in overpredicting the redshift of this galaxy {z^pec = 0.689) by > 4(T 
(-s^phot — S.BqIs, 3.4Q;g, and 3.4q;3, respectively). We point out that the large offset between the 
spectroscopic and photometric redshifts is possibly the result of source confusion, as two ra- 
dio sources (both with Spitzer and optical counterparts) fall near the center of the Bolocam, 
SCUBA, MAMBO, and SHARC 11 error circles, within 4" of each other. (Egami et al. 2004 
refer to the northwest and southeast radio sources as LE850.8a and LE850.8b, respectively.) 
The northwest radio source is believed by Lehmann et al. (2001) to be the counterpart to the 
ROSAT X-ray emission, who find a redshift of 0.974 using optical Keck spectroscopy. Using 
XMM-Newton observations, however, Ivison et al. (2002) conclude that the X-ray source 
corresponds to the southeast radio source. Indeed, the linear fit of a combination of ARP 
220 and MRK 231 ULIRG models to the Spitzer points coincident with the southeast radio 
source yields a 100% warm (AGN dominated) component. It is near this radio position that 
Chapman et al. (2005) find a spectroscopic redshift of 0.689. In fact, the X-ray emission 
observed with both instruments appears to fall between these two radio sources. Optical 
R-band images from both Yun et al. (2005), and Ivison et al. (2005) show multiple optical 
counterparts at the southeast radio source position. Furthermore, the spectroscopic position 
quoted by Chapman et al. (2005) appears to coincide with an optical source ~ 2" to the 
south of the southeast radio source, a source detected with four-color SDSS photometry (in 
addition to Yun et al. 2005 and Ivison et al. 2005 R-band photometry) and cataloged as a 
low-redshift galaxy. We conclude that it is possible that the submillimeter detections may 
either be suffering from source confusion from two or more galaxies, or that the Chapman 
et al. (2005) redshift corresponds to a source other than that of the southeast radio detection. 
If the latter is true, then the consistent redshifts predicted by the Laurent et al. /Wiklind 
(2005, 2003) and Spitzer IR models may further point to the accuracy of these photometric 
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redshift techniques. 

6. Discussion 
6.1. IR Spectrum 

SEDs over the entire IR-radio spectral range of the five Bolocam galaxies with spectro- 
scopic redshifts are shown in Figure 10. These spectra have the same redshift, cosmological 
dimming, and normalization corrections as in Figure 4. While four of the five galaxies have 
closely correlated spectra in the submillimeter region of the spectrum, the infrared spectra 
{Spitzer 3.6, 4.5, 5.8, and 8.0 /xm IRAC and 24 fim MIPS observations) exhibit a large degree 
of dispersion. This dispersion may be the result of several things: 1) Because the spectra 
have been normalized to a T = 40 K {P — 1.6) spectrum at their 1.1 mm Bolocam fiux 
densities to account for intrinsic brightness variation between the galaxies, the normaliza- 
tion will result in an artificial reduction in the submillimeter fiux density dispersions. This 
effect is not likely to dominate, as the flux density normalization has a ~ 10% effect on 
the flux densities of the galaxies. 2) The Spitzer detection associated with Bolocam source 
14 lies systematically low compared to the other three galaxies well-modeled by a T=40 K 
dust spectrum. While these Spitzer observations (Egami et al. 2004) are from a different 
data set than the remaining Spitzer observations (this work), it is unlikely that their flux 
densities are systematically uncertain by nearly an order of magnitude. 3) The 40 K dust 
temperature model of the submillimeter portion of the spectrum for Bolocam sources 5, 8, 
14, and 16 assume a single dust temperature for each of the four sources. The temperature 
fit parameter is somewhat degenerate with other fit parameters, including the critical fre- 
quency, fo, where the optical depth of the dust is unity. Thus, if Bolocam source 14 has 
a lower characteristic dust temperature, then the infrared portion of the 40 K model will 
significantly overestimate the infrared flux density. Chapman et al. (2005) estimates the 
temperature of Bolocam source 14 to be 33 K from two photometric points (850 fim SCUBA 
and 1.4 GHz VLA radio observations) and the dust SED templates of Dale & Helou (2002). 
This temperature uncertainty hkely dominates our uncertainty in matching the infrared flux 
densities. 4) In addition to heating by the ultraviolet and optical flux density from young 
stars associated with ongoing star formation, the thermal dust emission responsible for the 
bright submillimeter flux densities may be contributed to by an energetic AGN. While not 
dominating the total bolometric output from the galaxy, they may have a non-negligible 
(20%) contribution (Alexander et al. 2005). If this is the case, then the shape of the infrared 
continuum may be vary according to the relative contribution of star formation rates for 
these galaxies. 5) Another possible explanation for the larger dispersion in Spitzer infrared 
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flux density as compared to our single dust temperature model may be from the fact that 
the infrared flux densities trace separate epochs of star formation within the galaxy. It is 
plausible that the dust heated by the ultraviolet and optical flux density from from current 
star formation is not well correlated to the current infrared flux density of older stars (from 
previous star formation). Furthermore, models of UV to millimeter emission of star clusters 
embedded in optically thick giant molecular clouds (GMCs) suggest that the near-infrared 
to far-infrared portion of starburst galaxy SEDs vary considerably with age of the starburst 
(Efstathiou, Rowan- Robinson, & Siebenmorgen 2000). 

6.2. Radio Spectrum / FIR-Radio Correlation 

The composite SED (including the radio points) of the 5 Bolocam sources with spectro- 
scopic redshifts is shown in Figure 10. Two interesting observations can be made about the 
radio continuum emission associated with these galaxies: (1) the radio continuum is lower 
(or the submillimeter continuum is higher) on average than the well established radio-FIR 
correlation for the local universe by Yun, Reddy, & Condon (2001, cf. § 5.3); and (2) hke the 
infrared flux densities, the 6 and 20 cm VLA radio flux densities for the four galaxies with 
closely correlated spectra in the submillimeter region of the spectrum show a large degree of 
dispersion (factor of 5 in 20 cm flux density). This scatter is much larger than the quoted 
uncertainties of the VLA radio observations (which constrain the 20 cm radio flux density 
of each galaxy to better than 20%) This degree of dispersion is also larger than the factor 
of 3 scatter in the radio-FIR correlation seen among the FIR selected galaxies in the local 
universe. In fact, while the spectroscopic redshifts are similar for the four galaxies (2.1 < z < 
3.0), varying the 20 cm flux density over the observed range causes the best fit photometric 
redshifts of Yun & CariUi (2002) to vary from z — 2.7 to 4.8. This dispersion undoubtably 
contributes to the large errors of the photometric redshifts discussed in § 5. 

Deep radio continuum imaging using the VLA is a technically challenging task, and the 
disparate 20 cm flux densities by ~ 50% reported for Bolocam sources 2 and 6 (cf. § 4.2) 
exemplifles the difficulty of the photometry at radio wavelengths. Most systematic noises in 
interferometry tend to suppress the brightness of astronomical sources, and part of the lower 
radio continuum flux density might be related to the imaging and photometry problems. 

It is important to note that due to the low S/N ratios at which the sources have been 
detected in the submillimeter wavebands (< for the vast majority of SHARC, SCUBA, 
Bolocam, and MAMBO detections), flux bias (Laurent et al. 2005) plays a major role in 
overestimating the flux densities at these wavelengths. While this factor indeed results in a 
systematic shift of the entire submillimeter portion of the spectrum to higher flux densities. 
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it is unlikely that the magnitude of this effect (~ 20% for the Bolocam flux densities) could 
fully account for the lower than expected (from the FIR-radio correlation) radio continuum 
flux densities. Calibration errors may also contribute to a systematic overestimate of the 
submillimeter flux densities. Finally, as two or more independent radio sources are found 
within the Bolocam error circle in 8 out of 17 cases in Figure 2, source confusion or source 
blending may also contribute to the apparently lower radio continuum flux if only one radio 
source is identified as the counterpart. 

The lower radio continuum flux density and the larger scatter may also reflect an actual 
breakdown in the radio-FlR correlation. Invcrsc-Compton losses for the high energy cosmic 
rays responsible for the synchrotron radiation is thought to be significant at z > 2, and a 
possible breakdown in the radio-FIR correlation has been considered previously (see Condon 
1992, Carilli & Yun 2000). This effect is demonstrated in Figure 7, in which models of local 
ULIRGs Arp 220 and Mrk 231 systematically overestimate the radio flux densities of these 
submillimeter galaxies. Higher quality data on a larger sample of high redshift systems are 
needed to examine the importance of inverse- Compton loss and the possible breakdown in 
the radio-FIR correlation. 

X-ray heating of the circum- nuclear gas and dust is an important source of luminosity in 
the far-IR if a luminous AGN is present (Maloncy, Hollcnbach, & Tielens 1996). Radio-quiet 
AGN FIR emission could reduce the 1.4 GHz flux density with respect to the FIR heating. 
Alexander et al. (2005, and references therein) make the case using X-ray detections and 
spectral indices that many, perhaps most, submillimeter galaxies have AGN but that they 
are not bolometrically important. However, the statistics of X-ray detected submillimeter 
galaxies for which hard/soft ratios can be measured is not large and it cannot be ruled out 
that many submillimeter galaxies are Compton thick (Nh > 1.5 x 10^^ cm~^). We con- 
clude that the generally low radio flux densities in our sample could be due to small number 
statistics, source confusion, or generally depressed radio emission, perhaps due to quenching 
of high energy cosmic rays, although radio-quiet, Compton-thick AGN contributions to the 
dust heating cannot be ruled out. There are new AGN versus star formation spectral diag- 
nostics emerging (Ivison et al., Egami et al.), and it is possible that to deflnitively settle the 
issue may ultimately require ALMA, Constellation-X, and interferometric FIR spectral line 
diagnostic capability. 

6.3. Inferred Luminosities 

To obtain the intrinsic bolometric luminosities of the flve Bolocam galaxies with spec- 
troscopic redshifts, the spectra of § 4.4 (based on the Laurent et al. 2005 model) were 
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integrated for each galaxy in their respective rest frames. The resulting bolometric lumi- 
nosities are listed in Figure 4. The four galaxies well-modeled by a 40 K dust spectrum have 
luminosities ranging from L = (1.0 — 1.2) x 10^"^ Lq. The lower redshift galaxy (Bolocam 
source 17) has an inferred luminosity two orders of magnitude lower (L = 1.3 x 10^^ Lq). 
If the spectroscopic redshift of 0.689 does not apply to this galaxy and it instead lies at z 
— 3.4 (the photometric redshift predicted by the Laurent et al. 2005 / Wiklind 2003 and 
Spitzer IR models), then its luminosity of L = 8.2 x 10^^ Lq agrees well with the others. 

6.4. Stelleir and Dust Mcisses Implied from the Integrated Submillimeter 

Luminosities 

Three hundred and fifty micron observations combined with far-infrared and millimeter- 
wavelength photometry enables accurate measurements of luminosity for galaxies near z = 2 
because no interpolation across the peak of the SED is required. Characteristically, integra- 
tion of the SEDs of the galaxies in our sample from far-infrared to millimeter-wavelengths 
yields luminosities of ~ 1 x 10^^ Lq. Assuming: 1) that the luminosity derives from star 
formation (young stars, which may overestimate stellar masses due to contribution from in- 
termediate mass giants), 2) a characteristic (Salpeter) form of the initial mass function (IMF) 
from Chabrier (2003), and 3) all of the optical and ultraviolet radiation is reprocessed to long 
wavelengths by dust, enables the stellar mass content of the galaxies to be approximately 
estimated. We adopt the M^'^ luminosity function of Demircan & Kahraman (1991). Be- 
cause the luminosity function is so steep, the derived mass depends strongly on the assumed 
upper mass limit of integration. The lower limit of integration is not well constrained by the 
data, although masses less than 0.3 Mq are not hkely to dominate the mass because the IMF 
flattens considerably at low mass. Lower mass limits of 0.7 and 1.0 Mq could be relevant 
because: 1) Dwek et al. (1998) argued that a Salpeter IMF for z > 1 galaxies cannot extend 
much lower than this without producing too many low mass stars that would be present 
today and 2) Chabrier (2003) suggests (with caution) that the high-z IMF could cut off > 1 
Mq based on multiple circumstantial lines of evidence. Varying mi from 0.3 to 1.0 Mq and 
limiting ttIu to < 50 Mq yields a mininum stellar mass of 10^" Mq and a maximum of a few 
X 10^^ Mq, consistent with the stellar mass content of large elliptical galaxies, as previously 
pointed out by many authors (e.g.. Small et al. 2002; Lilly et al. 1999). This range must be 
considered an upper limit because AGN could be responsible for some of the dust heating. 

For all of the galaxies with secure 350 /xm detections, especially those with Chapman 
et al. (2005) spectroscopic redshifts, it is clear that the Bolocam 1.1 mm observations lie on 
the optically thin Ray leigh- Jeans side of the SED, and therefore enable dust mass estimates. 



The flux density, S^,^ of a galaxy at an observed frequency, is related to the dust mass, M, 

by 



where Di is the luminosity distance to redshift z, Bi,/{T) is the Plank function evaluated at 
the emitted frequency, z/', and Ki, is the dust opacity. Using the range of observed Bolocam 
flux densities (4.0 - 6.8 mJy), assuming a redshift of 2.1, and applying the dust cross section 
of K,,^ — 12 A cm^/g (Ossenkopf & Henning 1994) most relevant for high mass star formation 
(high gas density and thin ice mantle model), leads to dust masses of 3 - 5 x 10^ Mq (2-3 
X 10^ Mq) for a dust temperature of 40 K (50 K). Blindly applying a dust-to-gas mass ratio 
of 1/100 implies gas masses of 3 - 5 x 10^° M©. These gas masses arc comparable to those 
from a recent sample of 8 submillimeter galaxies of Gcnzcl et al. (2004) and Neri et al. (2003) 
which yield median molecular gas masses (from CO emission) of 2.2 x 10^*^ Mq and 2.8 x 
10^° Mq, with median dynamical masses of 1.1 x 10^^ Mq and 6.2 x 10"'^'^ Mq (assuming the 
most probable inclination angle of sin i = 2/n), respectively. These gas mass estimates are 
uncertain to at least a factor of a few due to 1) the Bolocam flux density bias, which causes 
the measured flux densities to be overestimated by 10 - 30%, 2) our assumed values of k, 
and T, which may vary by a factor of a few and ± 20 K, respectively, and 3) increasing our 
assumed redshift of 2.1 (the mean of the 5 Chapman spectroscopic redshifts of the Bolocam 
galaxies) to z = 2.4 increases our calculated gas mass by 30%. Nevertheless, taking these 
factor into account still imply that a considerable fraction of the mass could already by in 
stars and substantial gas remains for star formation. This major epoch of galaxy formation 
at approximately z = 2 is consistent with the conclusions of Fontana et al. (2004) from 
spectral fltting of a sample of 500 elliptical galaxies at 0.2 < z < 2.5 that approximately 
35% of eUiptical galaxy stellar mass was assembled hy z — 2 and approximately 80% by z 
= 1. 



7. Conclusions 

We have obtained 350 /xm SHARC II observations toward galaxy candidates from the 
Bolocam Lockman Hole survey. The Lockman Hole has rich, deep, multiwavelength obser- 
vations enabling detailed studies of galaxies. The 350 //m photometry is near the peaks 
of the SEDs of galaxies with characteristic temperatures of ~ 50 K and redshifts of 2; ~ 
2 to 3. They therefore enable measurements of luminosities and estimates of temperatures 
and photometric redshifts without interpolating over the peak of the FIR thermal SEDs. 
Seven galaxies detected at 1.1 mm with Bolocam were detected at 350 /im, two of which 
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have two 350 //m counterparts; these were combined with two 350 jim detections from the 
survey of Kovacs et al. (2005), bringing the total number of Bolocam galaxies detected with 
SHARC II to nine. Two additional galaxies not associated with the Bolocam sources were 
also detected. The SHARC II detections range in significance from 3.0 a to 6.8 a, with flux 
densities ranging from 14 mJy to 64 mJy. 

Wc combined our observations with 850 /im and 1.2 mm photometry from the literature 
to fit the submillimeter/millimeter-wave spectra to thermal dust models. We found that two 
models with significantly different dust temperatures (40 K and 68 K) and spectral indices 
P (1.6 and 1.8, respectively) yielded similar quality fits owing to the degeneracy in T and 
/3, rendering them indistinguishable without better SED sampling. However, there is little 
consequence of the degeneracy to the derived luminosities, photometric redshifts, and dust 
masses within the statistical uncertainties. Five of the galaxies have spectroscopic redshifts 
in the literature, with four ranging from 2; = 2.1 to 3.0 and one a,t z — 0.689. The four high-z 
galaxies have luminosities of (1.0 - 1.2) x 10^^ L©, while the z = 0.689 galaxy is best fit by 
a, 20 K, P — 1.0, spectrum with a much lower luminosity: 1.3 x 10^^ Lq. (Given the source 
confusion in the optical and radio, along with consistent photometric redshifts, we suggest 
that the z = 0.689 spectroscopic redshift of Bolocam source 17 may be a misidentification.) 
The characteristic dust masses for the four high- 2; spectroscopic galaxies are 4 x 10^ Mq, 
implying gas masses of 4 x 10^° Mq. The dominant uncertainties in this estimation are the 
dust opacity and the gas-to-dust conversion factor, which make the estimation uncertain to 
a factor of a few. Assuming a Salpeter IMF and that the submillimeter emission derives 
completely from star formation yields stellar masses of 10^° to a few times 10^^ Mq, broadly 
consistent with the stellar content of modern-day elliptical galaxies. 

The photometric redshifts of the full sample of seven galaxies span the range of 2; = 2.0 
to 2; = 4.3, with statistical uncertainties of Az — 0.3 to 0.6 (1 a). Photometric redshifts 
utilizing composite radio/FIR spectra representative of local star- forming ULIRGs yields 
systematically higher redshifts, on the order of Az = 1. For the four galaxies with optical 
spectroscopic redshifts the anomolously high redshifts arise from systematically low 1.4 GHz 
observed fiux densities. The discrepancy could arise from small number statistics, inverse- 
Compton losses of high energy cosmic rays off the CMB, heating by radio-quiet AGN, or 
suppressed synchrotron emission from supernova remnants in the unusually luminous galax- 
ies. For comparison, photometric redshifts derived using only Spitzer MIPS and IRAC data 
points yielded slightly more precise and accurate redshifts than the submillimeter/millimeter- 
wave data alone, with discriminatory power between heating by AGN and star formation 
(albeit with limited bolometric luminosity constraints). 

We thank Attila Kovacs for providing us with SHARC II detections of Bolocam sources 
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Fig. 8. — Histogram of photometric rcdshift errors for various models compared to the 
spectroscopic redshifts of Chapman et al. (2005). Yun & Carilh (2002) assume the far- 
IR to radio correlation while both Wiklind (2003) and Laurent et al. (2005) model the 
submillimeter and millimeter portions of the spectrum. The lower right frame is a fit of the 
near- and mid- infrared points to models of Arp 220 and Mrk 231 (Egami et al. 2004). 
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Fig. 9. — Flux density ratio between various wavebands assuming the submillimeter galaxy 
model SED of Laurent et al. (2005). The model is based on the composite SED of 
nearby dusty IRAS galaxies, high-redshift submillimeter galaxies, gravitationally lensed high- 
redshift galaxies, and high-redshift AGNs (Blain et al. 2002, and references therein). The 
Laurent et al. (2005) model assumes dust parameters of T = 40 K, i/q — 3700 GHz, /3 = 
1.6, and a power law of oc to model the hotter dust components on the Wien side 
of the spectrum. The shaded line represents the ratio of the SHARC 350 fim and Bolocam 
LI mm wavebands, with the thickness corresponding to the difference between the Laurent 
et al. (2005) and Wikhnd (2003) model SEDs. 
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Fig. 10. — Composite SED for the five Bolocam galaxies with spectroscopic redshifts (Bolo- 
cam 5, 8, 14, 16, and 17). The observed SEDs have all been redshifted to z = 2.0 using 
the spectroscopic redshifts, corrected for cosmological dimming, and their Bolocam 1.1 mm 
fluxes normalized to the Laurent et al. (2005) model. The solid line represents the Laurent 
et al. (2005) model (T = 40 K, /5 = 1.6) based on the composite SED of Blain et al. (2002, 
and references therein). The dashed line represents the same model, with T = 20 K and (3 
= 1.0. 



A. Multiwavelength Coverage 



Observations at both longer and shorter wavelengths are necessary to characterize high 
redshift submillimeter galaxies. Radio identifications provide precise astrometry which al- 
low optical and infrared counterparts to be identified. The radio observations detect the 
synchrotron emission (presumably from high energy cosmic rays associated with supernova 
explosions, see review by Condon 1992), which is known to be correlated to the far- infrared 
luminosity (and is thus correlated to the Bolocam flux density, which traces the dust emission 
that is responsible for the infrared luminosity) • The Bolocam and radio points may then be 
used to obtain temperature and redshift estimates. Multiple submillimeter wavelengths are 
essential in confirming the Bolocam sources and allow for a photometric redshift measure- 
ment that is independent of the radio-to-far-IR correlation (i.e. using only the submillimeter 
portion of the spectrum for an assumed temperature). Furthermore, a third photometric 
redshift estimate using the 1.6 //m bump in the stellar emission spectrum is possible using 
the near- and mid-infrared portion of the spectrum. In this appendix we briefly discuss the 
multiwavelength coverage (submillimeter, radio, infrared, optical, and X-ray) of the Lock- 
man Hole and the coincident detections of each survey. Figure 2 presents an overview of the 
Bolocam detections, including all of the submillimeter, radio, and infrared counterparts. 



A.l. Submillimeter Surveys 

Four submillimeter / millimeter surveys of the Lockman Hole region are considered. 
Three of the surveys (Bolocam, SCUBA, MAMBO) represent blank-fleld surveys, while the 
SHARC 11 survey covers only the 17 Bolocam sources. The surveys include: 

1) Entire Bolocam Lockman Hole survey at 1.1 mm (Laurent et al. 2005). The survey 
covers 324 arcmin^ down to an RMS of 1.4 mJy/beam and includes 17 source candidates (> 
3 a) ranging in flux densities from 4.0 to 6.8 mJy. Extensive simulations and jackknife tests 
yielded 6 expected false detections. The Bolocam beam size is 31" and observations were 
taken without chopping while scanning at 60"/s. The Bolocam survey includes the entire 
regions surveyed by the pubhshed 8 mJy 850 fj,m JCMT SCUBA (Scott et al. 2002) and 
1.2 mm IRAM MAMBO (Greve et al. 2004) surveys. Using a model SED based on nearby, 
dusty, star-forming galaxies (see § 5) gives relative flux densities of 1 : 2.0 : 0.8 and relative 
RMSs of 1 : 0.9: 0.6-1.4 for the Bolocam, SCUBA, and MAMBO surveys, respectively, for 
a galaxy redshift oi z — 2.4 (with the range given for MAMBO due to nonuniform noise). 

2) JCMT SCUBA 8 mJy survey (Scott et al. 2002). The survey covers 122 arcmin^ (in 
the Lockman Hole) to an RMS of 2.5 mJy/beam with a 14" beam, and includes 36 source 
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candidates (> 3 a). Five sources of the Scott et al. (2002) catalog (LE850.9, 10, 11, 15, 20) 

were retracted by Ivison ct al. (2002) on the basis of large agso/im values (and lack of radio 
identifications). The SCUBA survey detects 6 of the 8 Bolocam detections in the overlap 
region between the surveys (Bolocam 5, 8, 13, 14, 16, and 17). 

3) 1.2 mm IRAM MAMBO survey (Greve et al. 2004). The survey covers 197 arcmin^ 
(in the Lockman Hole) to an RMS of 0.6-1.5 mJy/bcam, with a 10.7" beam. The survey 
includes 23 sources (> 3 a). The MAMBO survey detects 7 of the 11 Bolocam detections in 
the overlap region between the surveys (Bolocam 1, 4, 5, 8, 14, 16, and 17). 

4) 350 iim SHARC II observations of the Lockman Hole (this work - See § 2). The 
instrument beam is 9" and the observations were taken with a Lissajous scan pattern with 
imiform coverage regions of 95" x 18" to an RMS of 6-18 mJy/beam (with deeper integrations 
taken for fainter sources). A total of 14 sources were detected (3 by Kovacs et al. 2005), 11 
of which coincide with 9 Bolocam sources. (An additional 3 SHARC II sources were detected 
outside of the Bolocam positional error circles.) 

To summarize, the majority of the Bolocam sources within the regions covered by other 
submillimeter surveys were detected (SCUBA detects 6 of 8 Bolocam detections, MAMBO 
detects 7 of 11, and SHARC II detects 9 of 17). Seven of the Bolocam source candidates 
were not detected by any other submillimeter survey (Bolocam 7, 9, 10, 11, 12, 13, and 
15, each of which was detected at < 3.5 cr), which is consistent with the 6 false detections 
(Poisson distributed) expected from simulations. These results point to the robustness of 
the remaining Bolocam galaxy candidates. 

A. 2. Radio Detections 

Three independent VLA radio surveys were used to identify counterparts to the Bolocam 
sources. In addition, this work re-analyzes the Ivison et al. (2002) data set to obtain a 
complete source catalog over the entire Bolocam source catalog. The radio surveys include: 

Published Surveys: 

1) VLA radio source list at 1.4 GHz of the SCUBA 8 mJy survey sources (Ivison et al. 
2002). The Lockman East radio survey reached an average noise level of 4.8 fiJj beam^^ with 
a 1.3" FWHM circular beam. Coincident sources (and 5a upper limits) are given for each of 
the 21 brighter (> 3a) SCUBA sources, with 20 positive detections over 14 of the sources. 
Five of the 21 SCUBA sources (LE850.9, 10, 11, 15, 20) are ignored on the basis of large 
as50fj,m values (and lack of radio identifications). Of the 8 Bolocam detections covered by this 
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VLA survey, 5 (Bolocam 5, 8, 14, 16, 17) coincide with 7 Ivison et al. (2002) sources. Ivison 
et al. (2002) also report continuum data at 4.9 GHz (using the VLA in the C configuration), 
with a resulting noise level of ~ 11 /iJy beam~^. Positive detections exist for 8 of the SCUBA 
galaxies, 4 of which are coincident with Bolocam sources (Bolocam 5, 14, 16, 17). 

2) A source list from the 6 cm VLA survey of the Lockman Hole (Ciliegi et al. 2003). 
The survey has a noise level of ~ 11/^Jy and covers the northeast quadrant of the Bolocam 
survey, including 7 of the Bolocam sources. A total of 63 sources were detected at > 4.5 cr, 
with a radio counterpart for only one of the Bolocam sources (Bolocam 8). 

Unpublished Surveys: 

3) A source list of 20 cm VLA observations was obtained from Yun et al. (2005). The 
observations were broken into four fields, with the northeastern section being deeper (10 
/iJy/beam RMS) than the other three fields (15 //Jy/beam). The observations cover the 
entire Bolocam field, with a total of 640 sources > 3.5o". The source list contains counterparts 
for 10 of the Bolocam sources (Bolocam 1, 2, 3, 5, 6, 8, 9, 13, 14, 17). 

4) Re-analysis of the VLA map from Ivison et al. (2002) by Biggs & Ivison (2006). The 
survey reached a median noise level in the locality of the 17 Bolocam sources of 4.4 //Jy 
beam~^. Radio sources were extracted from the VLA map within 30" from the Bolocam 
source positions. Counterparts for 15 of the 17 Bolocam sources were detected (all except 
Bolocam sources 10 and 12). Seven of the detections are coincident with the Ivison et al. 
(2002) analysis (Bolocam 5, 8, 14, 16, and 17), with flux densities in good agreement (within 
the stated uncertainties). 

In total, 15 of the 17 Bolocam galaxy candidates are found to have at least one radio 
counterpart (all except Bolocam 10 and 12). As only the VLA observations from Yun et al. 
(2005) and Ciliegi et al. (2003) include an entire source catalog, we are limited to estimating 
accidental detection rates for these two surveys. The accidental detection rate is defined as 
the Poisson likelihood that one or more of these known radio sources, randomly distributed, 
fall within the 2 a positional confidence region of the Bolocam beam (see § 3) . The rates for 
the Yun et al. (2005) and Ciliegi et al. (2003) surveys are 17-24% and 8-11%, respectively. 
(The accidental radio detection rate is dominated by the Yun et al. 2005 survey due to its 
greater relative depth.) 
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A.3. Infrcired Detections 

Both published and unpubhshed Spitzer observations were used to obtain infrared coun- 
terparts to the Bolocam / radio sources. They include: 

1) Spitzer postage stamps and an extracted source list (this work) with the IRAC (3.6, 
4.5, 5.8, and 8.0 ^m) and MIPS (24 fim) instruments were obtained for 14 of the 17 Bolocam 
source candidates for areas within a 16" radius of the Bolocam beam centers. (Note that 
this radius is smaller than the 2 a positional errors of Bolocam.) Bolocam. LE. 1100. 2, 6, and 
12 fall outside of the region observed with Spitzer. Multiple Spitzer sources were detected 
within each of the 14 Bolocam positional error circles. 

2) ^pzfeer observations of SCUBA/VLA selected sources (Egami et al. 2004) with IRAC 
(3.6, 4.5, 5.8, 8.0 /^m) and MIPS (24 /i,m) in a 5' x 5' area of the Lockman Hole region. The 
MIPS 24 iim field of view is 5.4' x 5.4'. The total integration time was 300 seconds per pixel, 
with a 3 (J detection limit of 120 /xJy. The observed region contains the following Spitzer, 
VLA radio, SCUBA, and Bolocam correlations: 

• The field imaged by IRAC and MIPS contains 10 SCUBA sources of the Scott et al. 
(2002) catalog (LE850.1, 4, 7, 8, 10, 14, 18, 23, 24, 35). 

• Ivison et al. (2002) detected radio counterparts for 5 (LE850.1, 7, 8, 14, 18) of the 
7 secure (> 3.5 a) SCUBA sources at high significance and a radio counterpart for 
LE850.4 with lower significance. 

• Egami et al. (2004) finds an additional radio counterpart to LE850.4 upon re-examination 
of the VLA map. 

• Of these 7 SCUBA sources with a total of 9 radio components (LE850.1, 4, 7, 8a, 8b, 
14a, 14b, 18, 35), all are found to have 5piteer counterparts. 

• Spitzer counterparts are found only for the sources with radio detections. 

• Four of the 9 Spitzer detections (LE850.1, 8a, 8b, 8c) fall within 2 of the Bolocam 
positional error circles (Bolocam.LE.1100.14 and 17). 

A. 4. Optical Detections 

Optical counterparts to the Bolocam / radio sources were obtained from three deep 
surveys, including: 
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1) SDSS u, g, r, i, z photometry was obtained from the pubhc SDSS data release 3 (DR3) 
website^, including 7,148 objects in the Lockman Hole region. The 61 objects falling within 
the Bolocam positional error circles include 43 galaxies (extended) and 18 stars (point-like). 

2) A hst of objects with Johnson R-band (SUBARU) photometry was obtained from 
Yun et al. (2005) for each of the Yun et al. (2005) VLA 20 cm radio sources. The optical 
catalog consists of 1,031 objects down to a limiting magnitude of R = 24.5. 29 sources were 
found within the positional error circles of the 10 Bolocam sources covered by the catalog. 
As no radio sources were detected by Yun et al. (2005) for Bolocam sources 4, 7, 10, 11, 12, 
15, and 16, no optical coverage was available for these sources. 

3) Optical Johnson R-band (SUBARU) photometry was also obtained from Ivison ct al. 
(2005) for most of the Lockman Hole region (with no coverage for Bolocam sources 6, 9, and 
11). 50,297 sources were detected down to a limiting magnitude of R = 27. 950 sources lie 
within the positional error circles of the 14 Bolocam sources covered by the survey. 

A. 5. X-ray Detections 

Finally, two X-ray surveys were used to identify Bolocam sources with a possible AGN 
contribution. The X-ray surveys include: 

1) The ROSAT Ultra Deep Survey (Lehmann et al. 2001) source hst of the Lockman 
Hole was cross referenced with the positions of the 17 Bolocam galaxy candidates. The 
survey reached a flux density level of 1.2 x 10~^^ erg cm~^ s~^ in the 0.5-2.0 keV energy 
band, detecting 94 X-ray sources. Two of the sources coincide with the Bolocam positions, 
coincident with the eastern radio source of Bolocam 8 and in-bctwccn the two radio sources 
associated with Bolocam 17. Despite the positional offset, Lehmann et al. (2001) associate 
the X-ray emission coincident with Bolocam source 17 with the northwestern radio (and 
optical) source. An additional X-ray detection exists of SHARC H source 3 (in the field of 
Bolocam source 9) , although the SHARC H and X-ray positions lie well outside the positional 
error circle of the Bolocam source. 

2) XMM-Newton observations of the Lockman Hole (~ 100 ksec) were taken during the 
Performance Verification phase of the instrument (Mainieri et al. 2002), yielding 98 sources 
with more than 70 net counts (flux density limit of 1.6 x 10~^^ erg cm~^ s~^ in the 0.5-7 keV 
band). Cross referencing the source list with the 17 Bolocam sources yields a coincidence 
list identical to that of the ROSAT survey. 



^http://cas. sdss.org/dr3/eii/tools/search/rect. asp 
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3) A deep, (1 Msec) unpublished (Ivison et al. 2005) XMM-Newton survey of the major- 
ity of the Bolocam good coverage region (Bolocam 2, 6, 11, and 12 are not covered) confirms 
the above results with only one identification, coincident with Bolocam source 17. 



B. Previous Redshift Estimates of Bolocam Galeixies 

Several of the Bolocam galaxies with radio counterparts have had their redshifts deter- 
mined with either photometric or spectroscopic techniques. We briefly discuss the results 
of each of these redshift surveys in this section. The spectroscopic redshifts will form the 
basis of comparison to the photometric redshift techniques discussed in § 5. The previous 
photometric redshifts also prove useful in considering the merits of including 350 //m SHARC 
II observations when estimating redshifts of submillimeter galaxies. 

1) Chapman et al. (2005) obtained spectroscopic redshifts for a sample of 73 submil- 
hmeter galaxies, the largest set of spectroscopic redshifts to date. The sources are a subset 
of 150 SCUBA/ JCMT sources detected at 850 fim over seven separate fields. Requiring 
accurate positions for spectroscopic identification, 104 of the SCUBA galaxies have radio 
identifications from deep VLA radio maps at 1.4 GHz. Of these, a subset of 98 sources were 
observed with the Low Resolution Imaging Spectrograph (LRIS) on the Keck I telescope, 
resulting in the spectroscopic identification of 73 of these galaxies. Twelve of the sources 
correspond to the SCUBA sources in the Lockman Hole, 5 of which a coincident with Bolo- 
cam galaxy candidates (Bolocam 5, 8, 14, 16 and 17). Chapman et al. (2005) find redshifts 
of 2.611, 3.036, 2.148, 2.142, 0.686 for the five Bolocam galaxies, respectively. While spectra 
for each of the Bolocam galaxies are not included among the sample spectra. Chapman et al. 
(2005) comment that all of the above detections have multiple-line identifications and can 
therefore be considered robust. The subset of Chapman et al. (2005) redshifts with Bolocam 
detections {z = 2.1 ± 0.9) is consistent with the overall redshift distribution of the entire 
sample (with a median redshift of 2.2). 

2) Aretxaga et al. (2003) implement a Monte Carlo photometric redshift technique, using 
existing far-IR- radio multiwavelength data for 77 sources first identified at 850 fim (SCUBA) 
or 1.2 mm (MAMBO). The technique produces the redshift probability distribution for an 
individual galaxy by choosing an evolutionary model and generating a catalog of luminosities 
and redshifts based on template SEDs from local starbursts, ULIRGs, and AGN. Thirteen of 
the sources correspond to SCUBA sources in the Lockman Hole, including 5 Bolocam galaxies 
(Bolocam 5, 8, 14, 16 and 17 - the same subset with Chapman et al. 2005 spectroscopic 



Table 3. Summary of Multiwavelength Detections of Bolocam Galaxy Candidates 



Bolocam 


Bolocam 


SHARC II 


SCUBA 


MAMBO 


Yun 


Yun (deep) 


Biggs 


Ivison 


Ivison 


Cilicgi 


Spitzer 


Spitzer 


Spitzer 


Spitzer 


Spitzer 


SDSS* 


Yun° 


Ivison 


Notes 


Number 


1.1 mm 


350 /Ltm 


850 ixra 


1.2 mm 


20 cm 


20 cm 


20 cm 


20 cm 


G cm 


G cm 


3.6 lim 


4.5 M^n 


5.8 lim 


8.0 fim 


24 /jm 


u,g,r,i,z 


R 


R 






mJy 


mJy 


mJy 


mJy 






MJy 


MJy 


MJy 


MJy 


mJy 


mJy 


m.Ty 


m.Ty 


m.Ty 












6 8il 4 


n-!-Q ^ 




5 7il 




60itl3 


52it5 7 




X 




1(4) 


1(4) 


-[-) 


( ^ 










radio source 








X 








53i5.6 


X 


X 




1 


1 




1 






X 




N radio source 


2 


6 5zbl 4 


9n Q-l-'^ 9 
^u. y iizo. ^ 






'^7^-1- 




t)u / my .u 




X 


X 


X 


X 


X 




X 


1 0(0 1) 




1 ('7A\ 


IN radio source 
















90i9 7 




X 


X 


X 


X 


X 




X 




X 




S radio source 




6 Oil 4 


15 lis 8 








-y 


9>^u^o n 
zoy m M .u 




X 


X 


1(10) 


1(4) 




) 




1(0 2)*^ 


2<1 




NE SHARC II source 






14.0±3.5 
















^ 


1 


1 














vv oxirtJr\,i_/ 11 source 




5.2±1.4 






3.2±0.7 


" 




20±6.7 




X 




X" (14) X'= (13) 


X (1) 


X (1) 


X {6} 


1,0(3,2) 


^ 


1(41} 






1 -H Q 


35.6i9.7 


llzt2.6 


2 9i0 7*^ 






51zt9-3 




60db35 




1(10) 


1(6) 






1 (o\ 




3d 


1(47) 


E radio source 












X 


59itl3 


42=t6.5 


X 


X 




1 


1 


1 


1 


1 


1,0 


3 


1 


S radio source 










3. lihO.7 


X 




27i5.9 


X 


X 




X 


X 


X 


X 


X 




X 




W radio source 




,j 












29i7.8 


''^ 






X 


X 




^ 


^ 








N radio source 


6 


5.0±1.5 


^7.o±o.y 


X 


X 




X 


55i9.9 


X 


X 


X 


X 


X 


X 


X 


X 


1,0(4,0) 




X 


hj D-HARU 11, Hi radio source 


^[ 


II 




X 


X 


105i22 


X 


138ill 


X 


X 


X 


X 


X 


X 


X 


X 






X 


Hi bMARLj 11, W radio source 






DO.D±iO.O 


X 


X 


611i93 


X 




X 


X 


X 


X 


X 


X 


X 


X 


1.0 


1 


X 


W brlARO 11 source 




4.93:1.5 
















X 


X 


-{6} 


-(5) 


-(IJ 






-{-) 




(•7^ \ 




8 


4.8il.3 


34.5i9.3 


lU.y=tz.4 


4.8i0.6 






27di7.6 


29dill 






1(13) 


1(13) 


1 / ^ 


-\-) 








1(83) 


S radio source 












X 


46ihl2 


28it6.7 


24it9 






1 


1 




1 


1 




-(2) 


1 


N radio source 


„ 


„ 


„ 






X 


136±23 


157ib8.0 


X 


X 


58itll 


X 


X 


X 


X 


X 


0,1 


3 


1 


E radio source | 


9 


4.8±1.5 


_c 


X 


X 


60±18 


X 


68±11 


X 


X 


X 


1(13) 


1(14) 


l(-) 


l(-) 


1(5) 


-(I'.O) 


2 


X 


CO 


10 


4.7±1.4 










X 




X 


X 


X 


-(6) 


-(6) 


-(3) 


-(1) 


-(1) 


-(2,1) 


X 


-(60) 


11 


4.6±1.5 




X 


X 




X 


116±33 


X 


X 


X 


1(6) 


1(4) 


-(-) 


-C-) 


-C-) 


1,0(3,2) 


X 


X 


12 


4.6±1.4 


_c 


X 


X 




X 




X 


X 


X 


X 


X 


X 


X 


X 


-(3,0) 


X 


-(66)*= 


1 


13 


4.5±1.4 




_c 




63±18 


X 


63±7.2 


X 


X 


X 


-(7) 


-(7) 


-(1) 


-C-) 


-C2) 


1,0(3,0) 


3d 


1(58) 


NE radio source 














X 


27±7.9 


X 


X 


X 


1 


1 


1 




1 


0,1 


X 


1 


E radio source 














X 


30±7.6 


X 


X 


X 














X 




SE radio source 


14 


4.4±1.3 


24.1±6.0 


lo.sii.e*^ 


3.4±0.6 


76±19 


78±15 


73±7.3 


73±10 


56±37 




X= (7) 


X<= (3) 


x^^ (-) 


X^^ (-) 


X"^ (1) 


-(1,1) 




1(83) 


S radio source 














X 


34±7.2 


X 


X 




1 


1 


1 


1 


1 




X 


1 


W radio source 














46±12 


26±8.4 


X 


X 




X": 


X'= 


X'^ 


x"^ 


X"^ 




1 


1 


NE radio source 
















26±8.3 


X 


X 














1,0 


X 


1 


SW radio source 


15 


4.4±1.4 




X 






X 


50±7.7 


X 


X 




1(13) 


1(12) 


1(3) 


1(2) 


1(2) 


1,0(5,0) 


X 


1(60) 




16 


4.1±1.4 


44.0±9.0 


6.1±1.8 


2.8±0.5 


X 




49±8.5 


41±12 


32±22 




1(9) 


1(7) 


l(-) 


l(-) 


l(-) 


-(1,1) 


X 


1(112) 




17 


4.0±1.3 


15.5±3.9 


5.1±1.3 


2.4±0.6 






20±6.6 


22±11 






1(4) 


1(5) 


1(1) 


1(0) 


1(2) 


0,1(3,2) 


1 


1(91) 


NW radio source 














54±13 


B8±9.4 


58±12 


B7±32 




1 


1 


1 


1 


1 




3(4) 


1 


SE radio source 



*The two listed numbers correspond to the SDSS galaxy and star classifications, respectively. 

'^The Optical R band source list (Yun et al. 2005) includes only sources within 3" of the corresponding radio source. 

■^See source description in § 4.2 

'^^Tlic multiple Yun et al. (2005) R-band optical sources coincident with the radio source are close in both position (< l") and flux (< 0.2 mag), and are likely due to a lack of merging close sources in 

the catalog. 



Note. — An (X) corresponds to Bolocam sources not covered by the survey / source list. Items in parentheses correspond to the detection of sources within the Bolocam positional error circle, but not 
likely associated with the source. 
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redshifts). Five different evolutionary models were used, with the results hsted in Table 4. 
Given the stated uncertainties in the photometric redshifts, the Aretxaga et al. (2003) results 
are in general agreement with the spectroscopic measurements. Note that the photometric 
redshift for Bolocam source 17 is consistent with both the modified blackbody and near IR 
stellar bump models (see § 5.4), further supporting a possible misidentification by Chapman 
et al. (2005). 

3) The ROSAT Ultra Deep Survey (Lehmann et al. 2001) of the Lockman Hole reached 
a flux density level of 1.2 x 10^^^ erg cm^^ s^^ in the 0.5-2.0 keV energy band detecting 94 
X-ray sources. Spectroscopic identifications of 90% of the sources based on highly accurate 
positions (2" FWHM) with the High Resolution Imager (HRI) and Keck R-band images were 
performed with the LRIS spectrometer on the Keck telescopes, yielding a single spectroscopic 
redshift for the northwest optical component of Bolocam source 17. Multiple lines were 
detected (Mg H, [Ne V], and [O H]) giving a redshift of 0.974 (cf. § 5.4). 

4) Oyabu et al. (2005) spectroscopically identified 29 of the 44 far-infrared sources (with 
accurate 1.4 GHz VLA radio positions) detected in the 0.9 deg^ ISO survey of the Lockman 
Hole. Optical identifications of the radio sources were performed in both the I- and R- 
bands (using the University of Hawaii 88" and Subaru telescopes) , with optical spectroscopy 
performed with both the Keck II and WIYN telescopes. Five of the ISO sources fall in 
the region covered by Bolocam, with 2 sources coinciding with Bolocam sources (6 and 8). 
Oyabu ct al. (2005) find a redshift of 0.362 for the most westerly radio source of Bolocam 6 
and a redshift of 1.110 for eastern radio source of Bolocam 8. 

5) With Spitzer IRAC and MIPS, Egami et al. (2004) positively detect all 9 of the radio 
sources associated with 7 SCUBA galaxies, including 4 detections of two Bolocam sources 
(14 and 17). The photometric redshifts using only the near- and mid-infrared points are 2.6 
for Bolocam 14, and ~ 3 and 0.9 for the two radio sources near Bolocam 17. 



Table 4. Previous Spectroscopic and Photometric Redshifts of Bolocam Galaxies 
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